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ABSTRACT 

"The  eatimation  of  the  ratio  of  two  independent  normal  vari¬ 
ances  is  considered  under  scale  invariant  squared  error  loss 
function,  when  the  means  are  unkno%m.  The  best  invariant  esti¬ 
mator  is  shown  to  be  inadmissible.  Two  new  classes  of  improved 
estimators  are  obtained,  one  by  extending  Stein  (1964)  and  the 
other  by  extending  Brown  (1968).  Numerical  studies  are  presented 
to  indicate  the  percent  improvements  in  risk. 


INTRODUCTION 


Let  X..,  j  *  l,...,n.,  and  n.  >  6  be  random  samples  from 

ij’  ■’  ’  ’  i’  1  -  2 

independent  normal  distributions  with  mean  and  variance  o^, 

i  *  1,2.  We  consider  the  problem  of  estimation  of  the  variance 
2  2 

ratio  6  *  o,/o-.  In  fact,  our  discussion  allows  innediate  exten- 
12  “l  ®2 

sion  to  more  general  parametric  functions  like  '“x 

and  m^  are  arbitrary. 

This  problem  is  motivated  by  the  work  of  Stein  (1964)  and 

Brown  (1968).  Stein  (1964)  proved  that  for  a  single  sample  of 

2 

X^'s,  the  usual  estimator  is  inadmissible  for  estimating  o  under 
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the  aquered  error  loaa  by  proving  an  eatiaMtor  which  haa  aaialler 
riak  (expected  loaa).  Brown  conaidera  the  more  general  problem  of 
eatimating  a  acale  parameter  in  the  preaence  of  an  unknown 
location  parameter  under  bowl-ahaped  loaa  leading  to  a  different 
claaa  of  dominating  eatimatora.  In  thia  paper,  ve  extend  both  the 
Stein  and  Brown  argumenta  to  two  independent  normal  populations. 

We  use  the  scale  invariant  quadratic  loss  function  of  the  form 

L(e,6)  -  (0  -  6)^6 (1.1) 

Let  X.  -  nT  X...  S.  -  e'(X. .  -  X.)^and  T.  -  n.XT,  i  -  1,2. 

^  ^  ^  j-1  ^  111’ 

Then  (X^  ,X2 ,8^^  .S^)  is  a  version  of  the  complete  sufficient 

2  2 

statistic  for  (Cj • 

As  is  well  known  (Stein,  1964;  Brown,  1968),  the  best  scale 

“1  .  2 

invariant  estimator  (n.  1)  S.  is  inadmissible  for  o.  under 

11  1 

squared  error  loss.  In  Section  2,  we  pursue  further  the  ideas  of 

Stein  and  Brown.  Our  work  there  is  related  to  that  of  Brewster  and 

Zidek  (1974)  and  Strawderman  (1974).  Our  findings  enable  us  to 

develop  in  Section  3  several  estimators  of  6  which  dominate 

«0  ■  ,  (1.2) 

the  best  invariant  estimator  of  6  under  loss  (1.1).  A  brief 
presentation  of  the  results  of  Monte  Carlo  simulation  to  measure 
percent  improvement  in  risk  is  included. 

2  2 

In  concluding  this  section,  we  argue  that  if  '''  , 

i  ■  1,2  independent,  then 


6  - 

2.2 


(1.3) 


is  admissible  for  6  ■  (1.1),  provided  n^  ^  5.  This 


n . 


will  imply  that  if  C,  are  known  by  taking  U.  •  E^(X. .  -  C.), 

12  ^  ®  1  ij  1 


j-1 


2 


i  •  1,2,  the  pair  ^U^U^)  is  complete  and  sufficient  for 


>> 

»  % 
fir. 

L 

r 


and  6  it  an  admitsible  eatimator  of  the  variance  ratio.  It  alto 

iiBpliet  that  6^  in  (1.2)  it  an  adniaaible  eatinator  of  6  in  the 

data  of  rulet  baaed  upon  (Sj^,S2). 

The  admiaaibility  of  (1.3)  nay  be  argued  a traight forwardly 

from  Brown  and  Fox  (1974,  pp.  808**810).  In  particular,  letting 

W  ■  logdJj/Uj),  V  ■  Uj,  n  ■  log  6,  ^  joint  denaity 

of  U  and  V  in  their  form  (3),  p.  809.  Uaing  the  generalised  prior 

d^/^  producea,  by  elementary  calculation,  (1.3)  as  an  invariant 

Bayes  procedure.  Their  regularity  conditions  for  admissibility 

(a)  -  (d),  p.  808,  will  be  satisfied  if  n^  ^  5.  (Condition  (d), 

the  most  difficult  to  verify,  holds  by  appealing  to  Bro%m  (1966, 

Theorem  2.3.3,  p.  1105).)  We  remark  that  if  we  use  a  gamma  prior 

6  6“1 

on  <(i,  i.e.,  f  „(♦)  -  e  ’^/r(6),  the  sai&e  argument  will 

a ,  p  ^ 

produce 


nj  +  26  -  4  Uj  +  a 

'~^2  07- 


(1.4) 


as  an  admissible  estimator  again  if  n^  ^  5.  Thus  (1.3) 

is  seen  to  be  the  limit  of  admissible  Bayes  as  well  as  generalized 

Baye  s . 

2.  IMPROVED  ESTIMATORS  OF  POWERS  OF  VARIANCE 

In  this  section  we  study  the  problem  of  estimating  arbitrary 
powers  of  variance  of  a  normal  distribution  under  the  scale  invari¬ 
ant  squared  error  loss.  Suppose  that  X^,...,X  is  a  random  sample 
2  —  2  ^ 

from  N(C,o  ).  Although  £(X.  -  X)  /(n  1)  is  the  best  estimator  of 

2  —  2^ 

a  in  the  class  cI(X^  -  X)  under  scale  invariant  squared  error 
loss,  Stein  (1964)  sho%rad  that  for  any  fixed 


/e(X.  -  E(X.  -  X)^\ 

■  -ttt- j 


dominates  l(X.  -  X)  /(n^-l)  under  the  loss 

T/  2  .  ,2  .2-4 

L(o  ,a)  *  (a  -  a)  o 


(2.1) 


(2.2) 


In  what  follows,  for  convenience  ve  take  Cq  *  0  and  let 

S  -  I(X.  -  X)^,  T  -  nX^. 

1 

In  eatinating  aasuming  ECS^™)  exiats,  c  .  s'"  where 

0~i  gID 


n-l  ,iB 


r(H^  ♦  ■)2' 

r(2^  ♦  2m) 


la  beat  in  the  claaa  cS  under  the  loss 


, ^  2m  «  /  20  V  2  -4m 

L(o  ,a)  ■  (o  -  a)  o 


(2.3) 


Again,  c  ,  s'"  is  inadmissible  under  the  loss  (2.3).  Follo%n.ng 

n—  1  ,m 

Stein's  approach,  we  consider  a  class  of  estimators  of  the  form 

A 

3(S/(S+T))(S+T)'"  for  o  '"  where  6(*)  is  a  function  from  [0,1] 
to,*).  Using  loss  (2.3),  the  risk  of  such  estimators  depends 

2  2.  5 

only  upon  C  /o  (so  that  without  loss  of  generality,  we  set  =  1) 

and  may  be  written  (e.g.,  Strawderman  (1974),  p.  191) 


-  2  r(5  ♦L+2m)2^'" 

1  +  E  -  V2L,m^  ”7(iTT) 


r(j  +1+0)  1 

r(|  +1+20)  2'^ 


(2.4) 


where  the  expectation  is  over  the  joint  distribution  of  S/(S'*'T)  and 
L,  with  S/(S+T)  having  a  beta  distribution  Be(^^,  — ^ — )  and  L 


2  2 

having  a  Poisson  distribution  with  parameter  n^  /2o  .  Clearly, 

c  increases  in  L  for  m  <  0  and  decreases  in  L  for  m  >  0. 

n+2L,m 


Hence,  if  6(S/(S+T))(S+T) 


n— 1  ,m 


one  can  choose  for  m  >  0,  B*  *  min(B,  c  )  and  for  m  <  0, 

n  giD 

B*  *  max(B,  c  ).  From  (2.4)  we  see  that  B*(S/(S+T)) (S+T)'" 
n  giD 

dominates  c  ,  s'"  under  loss  (2.3).  We  state  this  result  as 

n-l  ,m 

2 

Theorem  2.1.  For  X,  ,...,X  a  random  sample  fromN(C,o  ),  in 
■"  '  ■■  —  1  n 

estimating  o®,  m  >  -  (n-l)/4  (to  insure  ^  finite)  under  scale 

invariant  loss  (2.3)  if 


(2.5) 


0\ 

:  S  )  dominates 

o® 

c  ,  S 

n-l  ,m 

n-l  ,m 

,  s"*)  dominates 

c  ,  S 

n-l  ,m 

n-l  ,m 

In  particular,  when  a  ■  -1,  if  n  ^  5,  Bax(^“,  2^)  dominate • 


(n-5)/S. 

Brown  (1968)  has  created  a  different  type  of  estimator  for 


0^  (i.e.,  m  ■  1)  to  dominate  S/(n'»'l).  In  the  case  of  normality,  his 


estimator  takes  the  form 


S/(n+l)  if  T/S  >  c 
a*S  if  T/S  <  c 


(2.6) 


where  a*  <  l/(n-*'l)  and  depends  on  c.  In  what  follows,  we  give  a 
convenient  expression  for  a*  and  relate  the  Brown  and  Stein 
approaches.  Brewster  and  Zidek  (1974,  p.  22)  look  more  generally 
into  the  relationship  between  the  Stein  and  Brown  approaches  in 
the  context  of  dominance  of  equivariant  estimators. 


Using  (2.4)  at  m  *  1,  any  estimator  of  the  form 
6(S,T)  -  6(S/(S+T))(S+T) 


(2.7) 


under  the  invariant  lose  (2.3)  has  risk  of  the  form 


R(a,6(S,T))  -  El6(^)  "  (n+2L+2)'^)‘^(n+2L)(n+2L+2)  ♦  p(C  /a^), 

(2.8) 

tdtere  p  ■  E2(n  +  2L  +  2) 

By  noting  that 


jl-'  d  .  (1  .  c)-* 


and  defining 


®d,a*^y^  “  ^o,d)^y^y  *  **^(d,-)^y^y 


where  I  is  the  indicator  function,  the  Brown  estimator  is 


(2.9) 


i.e.,  of  the  form  (2.7)  with  risk  as  in  (2.8). 


The  beet  choice  of  •*  in  (2.9)  nay  be  obtained  from  expreaaion 
(6.2)  of  Bro%ni  (1968).  After  aome  manipulation,  we  may  write  it 
explicitly  in  the  form 


a*  -  (n+1)'^ 


(2.10) 


where  I^(a,b)  is  the  incomplete  beta  function  (tee,  e.g., 

Abramowitz  and  Stegun,  1965).  In  fact,  for  our  case,  we  can  simplify 
the  argument  in  Theorem  4.1  of  Bro%m  by  considering  (2.8)  in  the 
space  of  (S/(S+T),L).  We  omit  the  details. 

We  note  that  the  Stein  improved  estiaiator  (with  ■  0)  may 
be  written  in  a  form  similar  to  (2.6),  i.e., 

.-1 


S/(n+l)  if  T/S  >  (n+1) 

(S+T)/(n+2)  if  T/S  <  (n+1) 
.-1  . 


(2.11) 


-1 


In  fact,  if  we  use  c  •=  (n+1)  in  (2.6),  then  d  *  (n+l)/(n+2). 
Looking  at  (2.10)  we  see  that  if  d  ie  much  smaller  than  (n+l)/(n+2) 

(so  that  c  is  larger  than  (n+1)  *),  a*  is  essentially  (n+1)”^  and 
(2.6)  is  essentially  S/(n+l).  But  also  if  d  is  larger  than  (n+l)/(n+2), 

i.e.,  nearly  1,  then  c  is  very  small,  so  that  we  will  take  a*S 
with  very  small  probability  and  again  (2.6)  is  essentially  S/(n+l). 
Hence,  c  in  the  vicinity  of  (n+1)  ^  seems  best  and  our  numerical 
studies  support  this. 

Brown  (196§)  discusses  improved  estimation  of  o’"  for  m  >  0. 

In  fact,  analogous  to  Theorem  2.1,  %ie  can  show 

2 

Theorem  2.2.  Suppose  X^,...,X^  is  a  random  sample  from  N(C,o  ), 
in  estimating  o  ,  m  >  -(n-l)/4  under  loss  (2.3) 

s“  if  T/S  >  c 


n-1  ,m 


(2.12) 


a  S 
m 


>m 


if  T/s  <  c 

dominates  c  ,  s’"  in  terms  of  risk,  where 

n~i,m 


•  V 


c  ,  1,  .(1/2, 

n-l,«n  1-d 


n+2m-l 


>/I^_^(l/2,  (2.13) 


Proof.  The  proof  is  essentially  that  of  Brown's  Theorem  4.1.  When 
m  <  0  the  inequality  analogous  to  his  expression  (4.18) 

It 

is  reversed.  Notationally  we  have  suppressed  the  dependence  of  a 

in 

upon  n  and  d. 

4r  * 

Remark  2.1.  lfm>0,a  <c  ,  .ifffl<0.a  >c  ,  This 
— - ■  m  n-l,m  m  n-l,m 

is  analogous  to  using  a  minimum  or  maximum  according  to  m  >  0  or 
m  <  0  as  in  (2.5).  Alternatively,  the  estimators  in  (2.5)  may  be 
written  in  the  same  form  as  (2.12),  i.e., 

c  ,  s“  if  T/S  >  c  , 
n-l,m  n-l,m 

c  (S+T)“  if  T/S  <  c  , 
n ,  m  n~ i ,  m 

Therefore,  the  discussion  following  (2.11)  suggests  taking  c  in 

(2.12)  in  the  vicinity  of  c  , 

n*x  iin 

Remark  2.2.  Bro%m  (1968)  has  considered  more  generally  the 
estimator  of  a™,  m  >  0,  in  distributional  families  when  o  is  a 
scale  parameter  in  the  presence  of  an  unknown  location  parameter 
under  bowl'shaped  loss. 


In  this  section  we  will  show  that  when  the  means  are  unknown, 

2  2 

the  best  invariant  estimator  of  the  variance  ratio  6  = 
inadmissible.  We  will  obtain  several  improved  estimators  of  6. 

From  the  notation  developed  in  Section  1,  it  follows  that  the 
best  invariant  estimator  of  6  (when  the  means  and  unknown) 


(n^-S)  Sj 

"i**  h 


(3.1) 


Theorems  3.1  and  3.2  show  that  the  obvious  estimators  created  either 

2  -2 

by  improving  upon  the  best  invariant  estimator  of  Oj^  or  of 
dominates  (3.1). 

Theorem  3.1.  Under  the  lose  (1.1),  6^  is  inadmissible  for 
6  and  is  dominated  by 


'AV ' 


■"•?*  'V*’  ^ 


*1  *  "‘"‘‘o-  ^ 


n.-5  S,+T, 


(3.2) 


*2  '  ““<‘o’  ITh 


O’  n^+1  Sj+T^ 


(3.3) 


Proof .  The  proof  is  essentially  that  of  Theoreia  2.1  using, 


for  example,  in  the  first  case  the  fact  that  follows  up  to  a 


constant  a  central  chi-square  distribution  independent  of  T^. 


Remark  3.1.  Theorem  3.1  could  clearly  be  extended  to  provide 
dominating  estimators  for  ^  omit  the  details  here 


Taking  m^^  =  m^  *  1,  we  note  that  for  estimators  of  the  form 
6l(Si/(Si+Ti))  62(S2/(S2+T2))(Sj^+Tj^)/(S2+T2)  the  risk  depends  only 
2  2  2  2 

upon  and  without  loss  of  generality,  we  may  take 


*  1  and  hence  6*1.  Observing  that 


S  *'2”^ 

^1  “  — g —  min((Sj^+Tj)/(nj+2)  ,Sj/(nj^+l)) 


we  see  that  Stein's  result  is 


S,+T,  S.  2  S, 

E(.in(^  ,  .  1)  <  E(^ 


1)^  VC, 


while  we  have  added 


S,+T-  S,  n,-5  -  S,  n  -5 

47  -  i  47-  «!■ 


Thus  have  the  following  remark. 


Remark  3.2.  More  generally  using  the  same  argument  as  in 
Theorem  3.1  assuming  the  expectations  exist. 


S  S 

E{6*(^-^)(S^+Ti)  -  1}^  <  E{B(^-4^)(Sj+T^)  -  1}^  VC^ 


if  and  only  if 

S  S 

E{B*(^-^)(Si+Ti)h(S2)  -1}^<E{6(~t]^^^Si*'''i^^^S2^  ''^1  ’ 


Similarly 

E(6*(  J  )(Sj»Ij)  *  - 

11'  i  ^‘6‘s  J  V^2> 

2  2 

2  2 

if  and  only  if 

S,  h(£LXS  ♦T  )'^ 

2  S.  h(S,)(S-+T,) 

^S.+T.^  n.+l  ■ 

1)  i=<6(s  .T  >  „.l 

Remark  3.3.  It  is  noteworthy  that  the  argument  for  Theorem  3.1 
and  Remark  3.2  depends  entirely  upon  using  S.  with  estimates  of  the 
Si 

form  j+Tj)  (or  vice  versa)  rather  than  using  S2  ♦  T2. 

Thus  we  cannot  show,  for  instance,  that  the  appealing  estimator 
_  S,+T,  S,  n*~4  n-“5 


S,+T,  S.  n*-4 

*6  ■  "‘"<^75-  •  •  -s7> 


(3.4) 


dominates  6^^,  6^  or  even  6^.  We  will  return  to  this  point  later. 

With  regard  to  the  Brown  estimators,  we  have 
Theorem  3.2.  Under  the  loss  (1.1),  <5^  is  inadmissible  for  6  and  is 
dominated  by 


if  Tj/S^  >  c 
if  T^/Sj  <  c 


(3.5) 


If  •  ^ 

60  If  T^/S^  >  c 

6  *  *  1 
2  ^  (^)  if  T,/S,  <  c 

n.+i  Sj  22 

where  a*  is  given  in  (2.10)  and  a**  E  a  ^  in  (2.13). 
Proof .  Consider  6^^.  We  may  write 


(3.6) 


and  appeal  to  Remark  3.2.  The  argument  is  similar  for  62. 

Remark  3.4.  Extension  of  Theorem  3.2  in  the  direction  of  Remarks 
2.2  and  3.1  is  straightforward.  We  omit  the  details. 


The  estimatort  5T  and  6.,  i  ■  1,2.  are  not  admiasible.  We 

next  provide  explicit  eatimators  which  dominate  them.  In  the 
process  we  clarify  the  problem  of  estimation  of  6.  Suppose  that 

“l  “  Sj/(nj+l),  a®  •  (Sj+Tj^)/(nj+2),  o®  «  a*Sj 
and  define  regions 

Aj  -  {Tj/Sj  >  (n^+l)“^  T^/Sj  >  (n^-Sr^} 

Aj  *  {Tj/Sj  <  (n^+1)"^,  T^/S^  >  (n^-S)"^  ^3 

A3  »  {Tj/Sj  >  (nj+l)"\  T^/S^  <  (n^-S)"^ 

\  =  {T^/S^  <  (n^+1)*^.  T^/S^  <  (n^-S)'^. 


Then  Figure  1  shows  the  previously  discussed  6's  and  several 


additional  ones.  In  looking  at  the  dominance  of  63  over  6^,  we 


can  show  that  is  always  accomplished  on  A2,  %*hereas  the  situation 


is  unclear  on  A, .  Similarly  the  dominance  of  6.  over  6.  is  always 


accomplished  on  A3,  again  with  the  situation  on  A^  unclear.  The 


following  theorem  gives  the  dominance  results  mentioned  after 
Remark  3.4. 

Theorem  3.3.  Under  the  loss  (1.1)  the  following  hold: 


S  SB  B 

63  dominates  63,  63  dominates  63, 


S  SB  B 

6^  dominates  6^  and  6^  dominates  6^. 


Proof .  Immediate  from  the  Lemma  1  of  the  appendix  and  the 
definition  of  the  6's. 


S  B 

Corollary  3.1.  6^  dominates  6^  and  6^  dominates  6^. 


Remark  3.5.  For  the  Brown  estimators,  one  can  use  general  C3, 


in  defining  the  A^'s  and  the  above  results  will  go  through. 


Remark  3.6.  The  argument  in  Lemma  1  of  the  appendix  fails  on  A. 


and,  in  fact,  there  is  no  best  choice  on  A,.  For  example,  the 

s  .  ,  ^  s  s 

intuitively  appealing  estimator  6,  in  Figure  1  suggests  on  A. 
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FIG.  1.  Estimators 

of  6  Defined 

by  Regions 

but  we 

g 

cannot  show  that  6.  dotainates  even  6^. 

o  0 

However  on  A, 

4 

( sup- 

pressing  superscripts), 

0261  <  0262  ^  “l®2'  suggests  that  a^Bj^  will  be  close  to  012^2’ 

that  is,  the  performance  of  6,  will  be  essentially  that  of  6c  which 

D  j 

does  dominate  6n'  Our  numerical  studies  show  that  6c  or  6,  is  near- 
u  j  b 

ly  always  best  supporting  a  middle  rather  than  an  extreme  choice  on 

Monte  Carlo  simulations  were  performed  to  obtain  risks  for  the 

estimators  in  Figure  1.  10,000  replications  were  used.  The  per- 

S  S  S 

centage  improvements  in  risk  with  respect  to  6q  for  6^,  6^,  6^, 

g 

given  in  Table  1  for  the  case  n^  =  n„  =  10, 
Dj*^*5*6  i>fc 

2  2  2  2 

=  1  and  a  range  of  Cj.  ^2’  expected  the  performance  of 

6^  is  close  to  that  of  6^  although  generally  a  bit  better.  The 


V-*-.  -v-.*- 
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percent  improvement ■  are  email  and  become  smaller  with  increasing  n^. 
Hotiever  they  are  greater  than  those  observed  by  Brown  (1968).  His 
brief  numerical  study  showed  for  a  single  variance  under  squared  error 
loss  a  maximum  improvement  of  1  to  2Z.  For  a  variance  ratio  we  are 
able  to  roughly  double  this.  Moreover,  the  simplicity  of  the 
Stein  estimators  encourages  their  use  particularly  for  small  n^.n^. 

TABLE  I 


Percentage  Improvements  in  Risks* over  6^ 


”1 

■"2  ■ 

10 

6® 

4 

A* 

(0.0) 

2.0 

4.0 

2.9 

4.1 

2.0 

3.5 

3.0 

3.5 

(0,.01) 

2.1 

4.0 

3.1 

4.0 

2.0 

3.8 

3.0 

3.7 

(0,.l) 

1.5 

3.2 

2.5 

2.9 

1.7 

3.5 

2.6 

3.3 

(0,1) 

0.3 

0.6 

0.4 

0.4 

0.5 

0.8 

0.6 

0.7 

(.01,0) 

2.5 

4.0 

3.2 

4.2 

2.4 

3.4 

3.2 

3.6 

(.01,  .01) 

2.2 

4.1 

3.1 

4.0 

2.0 

3.9 

3.0 

3.7 

(.01,.!) 

1.6 

3.4 

2.4 

3.1 

1.8 

3.5 

2.6 

3.3 

(.01,1) 

0.4 

0.7 

0.6 

0.6 

0.5 

0.9 

0.6 

0.7 

(.1,0) 

3.5 

4.1 

3.9 

4.3 

2.9 

3.3 

3.4 

3.5 

(.1,.01) 

3.3 

4.0 

3.8 

4.2 

2.8 

3.6 

3.4 

3.7 

(.l,.l) 

2.6 

3.0 

3.2 

3.3 

2.7 

3.3 

3.2 

3.4 

(.1,1) 

0.8 

1.0 

0.9 

0.9 

0.8 

1.0 

0.9 

0.9 

(1,0) 

3.1 

3.0 

3.1 

3.0 

2.4 

2.3 

2.3 

2.3 

(1,.01) 

2.9 

2.8 

2.8 

2.9 

2.7 

2.5 

2.6 

2.6 

(l,.l) 

2.2 

2.2 

2.2 

2.2 

2.4 

2.4 

2.4 

2.4 

(1.1) 

0.3 

0.3 

0.3 

0.3 

0.4 

0.4 

0.4 

0.4 

*The  largest  sample  standard  error  over  all  the  128  percentage 
improvements  was  less  than  .2. 


APPENDIX 

Lemma  1.  Under  the  notations  in  (3.7)  with  a,  either  a~  or  a?  and 
”  SB  2  2  2 

62  either  $2  or  621  the  following  inequalities  hold: 

(i)  El.  ...  (a  6  -1)2  <  El  (a  6-1)2  v5 

A2UA4  2  1  -  A2UA,  1  1  1 


El.  ^  El* 

A2  2  1  —  A2  i  A  A 


(A. 2) 
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Proof.  Simple  calculation  ahowa  that  (A.l)  ia  equivalent  to 

i  2“a2UA,®1<»1-»2>- 

Note  that  on  >  a^.  Thua  (A. 3)  becomea 

/  (Oj-  a2)daj^da2 

“l''“2  2 

■■■->■-= - 1  E(b:)/2E(6,). 

/  (a^-a2)dajda2 
“l"“2 

We  need  only  show  that 


(A. 3) 


E(Bp  ^  ECeJjTj/S^  >  (n^-S)"^ 
ECejT^/S^  >  (n^-S)"^ 


(A. 4) 


Define  the  density  pCs^)  of  $2  hy 


p(s2) 


^l^®2^^n2-l^®2^ 


Then  for  each  fixed  t2,  taking  expectations  w.r.t.  the  pdf  p(s2^* 
one  gets 

E  (6,)E  J  <  E  16-  I.-  f  ..  J.  (A. 5) 

pip  (0,in2-5)t2)  —pi  (0,(n2-5)t2) 

Taking  expectations  over  t2  of  both  sides  of  (A. 5)  and  simple 
manipulation  yields  (A. 4).  The  proof  of  (A. 2)  is  similar. 
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in  risk. 
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